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Abstract 

We complete the leading-log renormalization group scaling of the NRQCD La- 
grangian at 0(l/m 2 ). The next-to-next-to-leading-log renormalization group scaling 
of the potential NRQCD Lagrangian (as far as the singlet is concerned) is also obtained 
in the situation ma s 3> Aqcd- As a by-product, we obtain the heavy quarkonium 
spectrum with the same accuracy in the situation 111^ > Aqcd- When Aqcd <C ma 2 s , 

this is equivalent to obtain the whole set of 0(mai n+4: ^ ln n a s ) terms in the heavy 
quarkonium spectrum. The implications of our results in the non-perturbative situa- 
tion ma s ~ Aqcd are also mentioned. 
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1 Introduction 



Heavy quark-antiquark systems near threshold are characterized by the small relative veloc- 
ity v of the heavy quarks in their center of mass frame. This small parameter produces a 
hierarchy of widely separated scales: m (hard), mv (soft), mv 2 (ultrasoft), ... . The fac- 
torization between them is efficiently achieved by using effective field theories, where one 
can organize the calculation as various perturbative expansions on the ratio of the different 
scales effectively producing an expansion in v. The terms in these series get multiplied by 
parametrically large logs: In if, which can also be understood as the ratio of the different 
scales appearing in the physical system. Again, effective field theories are very efficient in 
the resummation of these large logs once a renormalization group (RG) analysis of them has 
been performed. This will be the aim of this paper for the cases of NRQCD |IJ and potential 
NRQCD (pNRQCD) [g §0. 

First, we will obtain the RG improved matching coefficients of the NRQCD Lagrangian 
at one loop and up to 0{l/m 2 ). Since, by construction, the matching coefficients of HQET 
are equal to the analogous ones of NRQCD, these can already be obtained from the literature 

0, ||. Therefore, only the four-heavy fermion matching coefficients need to be computed 
to obtain the complete leading-log (LL) RG improvement of the NRQCD Lagrangian at 
0(1 /m 2 ). We will perform such a calculation in this paper. For the spin-dependent four- 
heavy fermion matching coefficients there already exists a computation in Ref . || . We differ 
with their results. Our evaluation is relevant in the study of the situation ma s ~ Aqcd- 
On the one hand they could be used to improve the accuracy of phenomenological stud- 
ies or lattice simulations of NRQCD. On the other hand, this situation has been studied 
within an effective field theory framework in |7| where the matching between NRQCD and 
a Schrodinger like formulation has been achieved in a controlled fashion. In particular, it 
is possible to write the potentials as Wilson loops multiplied by the matching coefficients 
inherited from NRQCD. The obvious application is that the matching coefficients here com- 
puted are the ones that multiply the Wilson loops in the non-perturbative potentials. This 
is especially relevant now that the complete expression for the potential at 0(l/m 2 ) in terms 
of Wilson loops is available J7|, |[ [| . In particular, it would be welcome to have an updated 
evaluation of the lattice analysis of the heavy quarkonium spectrum made in Ref. |1 taking 
into account the complete 0(l/m 2 ) potential as well as the now known complete set of LL 
NRQCD matching coefficients. 

In the situation when ma s ^> Aqcd, the matching between NRQCD and pNRQCD, 

1. e. the computation of the potentials, can be done perturbatively. In this case ultrasoft 
gluons as well as the quark-antiquark in an octet configuration do exist at the matching 
scale between NRQCD and potential NRQCD producing further divergences. By taking 
into account these divergences as well as the divergences computed before we have obtained 
the next-to-next-to-leading-log (NNLL) RG improved pNRQCD Lagrangian (as far as the 
singlet is concerned). 

If we are in the situation Aq CD / (ma s ) 2 <C ma 2 , the leading solution of the spectrum cor- 
responds to a Coulomb-type bound state and the non-perturbative effects are corrections. In 

■"■We will use here the pole mass as the expansion parameter leaving aside in this paper any considerations 
about renormalons. 
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this situation, by using the previous result of the NNLL RG improved pNRQCD Lagrangian, 
we are able to obtain the heavy quarkonium spectrum with the same accuracy. 

If instead, we are in the situation Aq CD <C ma 2 , from our previous result, we are able to 
obtain the whole set of 0(ma| n+4 ' ln n a s ) terms in the heavy quarkonium spectrum. There 
already exists an evaluation |TI| within the vNRQCD framework |TT[ [T2|, |T^] of the RG 



improved Heavy Quarkonium mass when Aqcd *C ma 2 s . We agree for the spin-dependent 
terms (since we do for the spin-dependent potentials computed in Ref. |T2|) but differ for 
the spin-independent ones. 



2 NRQCD 

NRQCD has an ultraviolet (UV) cut-off vnr = {v p , v s } satisfying mv <C i>nr <C th. At this 
stage Up ~ u s . u p is the UV cut-off of the relative three- momentum of the heavy quark and 
antiquark. v s is the UV cut-off of the three-momentum of the gluons and light quarks. This 
does not seem to give problems at the order we are working at but one should be eventually 
careful upon the possible gauge dependence of this splitting. 

Indeed, the above cutoffs plus the matter content of the theory given below correspond to 



our definition of NRQCD. Within the threshold expansion framework [II] this corresponds to 
integrate out the hard modes of QCD in order to obtain NRQCD. Unfortunately, NRQCD 
already contains non-physical degrees of freedom for the phase space region it is aimed 
to describe, which implies that the terms in the Lagrangian do not have a unique size nor, 
therefore, power counting (to avoid this problem is one of the motivations for the construction 
of pNRQCD, which will be done in the next section). Nevertheless, this poses no problem 
for the NRQCD running considered in this section. 

The NRQCD Lagrangian including light fermions reads at 0(l/m 2 ) (up to field redefi- 
nitions) |, 1, § 

£ = C g + Ci + + C x + L,^ x , (1) 

where ip is the Pauli spinor that annihilates the fermion, \ is the Pauli spinor that creates 
the anti-fermion, iD = id$ — gA , zD = z'V + gA, 

£ g = -\G^ a G^ + d[^gf abc G%G b \G c ™, (2) 



g 2 g 2 

+ c 3 g-s E %7 M & qjJM + 4 %7 M 7s?i Qjl^qj , (3) 



J „ D 2 D 4 a B 

C, = ^{^ 0+ c k —+c i — + c F g— (4) 

(D-E-E-D) cr 

+CDg \-icsg 
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2 2 

+ c i^—2 E 1^oT a qi + ^/-j £ ^7 W ^75^ 

8m^ ~^ 8rrr i 

2 2 

8m/ 8?rr 
analogously for £ x and 

C+x = %V Wx + % WxW + %t T > x t T * x + %t T V^t T V X . (5) 

m Z m A m Z m Z 

We have also included the D 4 /m 3 term above since it will be necessary in the evaluation of the 
heavy quarkonium mass once the power counting is established. Moreover, we will consider 
that the kinetic term matching coefficients are protected by reparameterization invariance 
(cfe = C4 = 1) however, we will often keep them explicit for tracking purposes. 

The NRQCD matching coefficients are functions of u^r = {v p , v s }. Somewhat by def- 
inition, the matching coefficients of the bilinear in the heavy quark fields and of the pure 
gluonic terms are just functions of u s , i.e. c = c(u s , m) = c(u s ). In any case, it will explicitely 
come out from the calculation. The complete LL running of these matching coefficients in 
the above basis (@-fD have been calculated by Bauer and Manohar || in the (background) 
Feynman gauge! (some partial previous results already exist in the literature [|J). Therefore, 
in order to complete the RG running of the NRQCD Lagrangian we only need to compute 
the four-heavy-quark operators with LL accuracy. This will be our aim in the following. 

A procedure to get the v s dependence of the NRQCD matching coefficients is by using 
HQET-like rules in NRQCD (by this we mean to perform the perturbative expansion in 1/m 
prior to the computation of the Feynman integrals). In fact, in our case, at one loop, all the 
dependence of the matching coefficients is only due to v s since no v v dependence appears at 
one-loop, i.e. d{y v , v s , m) = d(v p , v s ) ~ d{y s ). This will be discussed below within pNRQCD. 

Formally, we can write the NRQCD Lagrangian as an expansion in 1/m in the following 
way: 

OO 1 

£NRQCD = £ ~^K°n7 ( 6 ) 
n=0 m 

where the above fields and parameters should be understood as bare and the renormalization 
group equations of the renormalized matching coefficients read 

„A a = Ba(a) . (7) 

The RG equations have a triangular structure (the standard structure one can see, for 
instance, in HQET RG evaluations, ie. for the Lagrangian ([I]) setting the heavy antiquark 
field to zero): 

v s -, — A = i?o(Ao), 
dv s 

v s - — Ai = £?i(Ao)Ai, 
dv s 



v s -^-^2 = -B 2 (2,i)(Ao)A 2 + -B 2 (i,2)(Ao)Ai . 



2 We thank C. Bauer for communication on this point. 
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where the different B's can be power-expanded in Ao (Ao corresponds to the marginal op- 
erators (renormalizable interactions)). For NRQCD we have Ao = a s and Ai = {ck,c F }, 

\2 = {c 9 1 ,c D ,c s ,{c a },{c hl },{d}}- 

At this stage, we would like to stress that we are working in a non-minimal basis of 
operators for the NRQCD Lagrangian. Consequently, the values of (some of) the matching 
coefficients are ambiguous (only some combinations with physical meaning are unambiguous). 
In particular, some of the matching coefficients could depend upon the gauge in which 
the calculation has been performed. Therefore, it is important to perform the matching 
calculation in the same gauge (at least for those operators which could suffer the ambiguity). 
We will further discuss this point latter on. 

The RG equations for the {c} in the Feynman gauge can be read from Bauer and Manohar 
results 0. Because of latter comparison, we only explicitely write the equation for cd, it 
reads 



d a s 

' dv* 47T 



AC A (2C A 32CA 2 10C A 2 8T F n f hl 

Cn — CZ C% H -C\ 

3 V 3 3 / k 3 F 3 1 



(9) 



where, T F = 1/2, Cf = N £ N 1 and C A = N c . The explicit expression for the Cn RG equation 
depends on the gauge!. 

The RG equations for the d's in the Feynman gauge are new and read! 



_d_ 

dv s 

v s —d vs = A(Cf - C A )a 2 s cl + -a 2 s C A c D , 
dv s 2 

d j Ca 22 fii\ 

Us dV = — Y asCp ' ^ ' 

These equations have been obtained by explicit computation in the Feynman gauge by Signer 
fl6|| within the threshold formalism |14| . We have obtained them by using the results of Ref . 



17J , which were performed in the Feynman gauge, plus doing the explicit calculation of 
the terms that depend on cr> in NRQCD in the Feynman gauge. This proves to be enough 
since the dependence on c| of Eq. flllD can be inferred from the results of Ref. [|17] once 
the dependence on is known (since the spin-dependent terms will depend on c F ). Both 



3 The renormalization group evolution of the one-heavy quark sector has also been done in a minimal 
basis in Ref. j^] by elliminating the operator multiplying c\ l . In that case the expression of co obtained in 
Ref. H is indeed gauge independent. 

4 For the record, we also display the non-equal mass case equations with the definitions d/m 2 — > d/{mim2}- 
The equations for d ss and d sv remain equal, for d vv one has to change c 2 F — > c^Cp^ and the equation for 
d vs reads 

d j ( A n 3C a\ 2 2 , 3 2 n (mi (2) m 2 5 2 2 /mi m 2 \ , . 

v s -7-d vs = 4(7/ — a s c k + -a s C A — c y D ' + — & D > - -c k a s C A + . (10) 
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calculations agree. Note that it was needed not to have v v dependence at one-loop in order 
the argument to go through. 

As we have mentioned we are not working in a minimal basis. This shows up in the 
ambiguity of the value of the matching coefficients of some operators. At the practical level, 
this means that they will depend on the specific basis of operators we have taken for the 
NRQCD Lagrangian and on the procedure used (in particular on the gauge). Therefore, if 
working in a non-minimal basis, one should be careful and do the matching using the same 
gauge for all the operators (or at least for those that are potentially ambiguous). 

For illustration, let us consider the case without light fermions. In this case, cd and 
d vs are ambiguous but not an specific combination (see Eq. (0)). In particular, cd could 
be absorbed by other matching coefficients by using a field redefinition ||. We can check 
these statements by doing the calculation in the Coulomb gauge. In this case we obtain (no 
non-trivial change is now required in the non-equal mass case for the RG equations of the 
four-heavy fermion matching coefficients) 



v s - — cd (Coulomb) = 

dVe 4:71 



3 -° D V 



u s — ^.(Coulomb) = [iC f -^) a y k . (12) 



One can see that, as far as the combination that appears in Eq. fl34|) is concerned, the 
physical result is unchanged. 

In the following we will use the Feynman gauge results for the NRQCD matching coeffi- 
cients. 

With the above results we have completed the RG equations of the NRQCD Lagrangian 
at one loop at 0(l/m 2 ). In order to solve these equations, we need the (tree- level) matching 
conditions of the matching coefficients at some matching scale. We choose as the matching 
scale m. The {c(m)} can be read, for instance, from ||. The tree-level matching conditions 
for the four-heavy fermion operators read 

d ss {m) = -3C f (c f - Tca s (m) , 
d sv (m) = C f (Cf - -£\ ■Ka s {m) , 
d vs (m) = 3 (c f - -y^ -na s {m) , 

d vv {m) = - (c } - iia s {m) . (13) 
We can then obtain the solution of the RG equations. We only explicitely display those 

ot s {m) . 



which are new or will be necessary later on (we define z = ^4^T ^ — ^ — V (^ n ) a s( iy s) hr( ^ ) 
c F (u s ) = z 



3- 

C A 

-C A 



c s (v s ) = 2z~^-l 
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cd\v s 



d s 
d„ 



9C A f + _ 2Ca C a + lQC f - 8T F n f 

9C A + 8T F n f \ AC A + AT F n f Z 2(C A -2T F n f ) 



+ 



-7C 2 A + 32C A C 



f 



AC A T F n f + 32C f T F n f 4Ti , raf /3 _ 2 cu/3 



A(C A + T F n f )(2T F n f -C A ) 



d s 
d r 



9C A 
(m) + AC f 



Z -2C A + 



20 32 C/ 
13 + 13 ~C~ A 



C f - 



[m] 



m 



G, 



57T 



7T 



/3o 



a s {m) 



00 



l-z 



- 1 



Po 



27C 2 , 



7T 



a a (m) 



9^ + 8^71/^ 

C A + 16C f - 8T F n f 



5C A + 4T F n f & 
~4C A + 4T F n f p -2C A 

2(C A -2T F n f ) V J 
-7C1 + 32C A C f - 4C A T F n f + 32C f T F n f 
4{C A + T F n f ){2T F n f -C A ) 



X 



3/3o 



3A, + <iT F rif - 2(7, 



/3 +4T F n / /3-2C A /3 



8T F nf 
' 9C A 



00 



d vv {v s ) = d vv (m) + 



Po ~ 2C A 

c A 



Po-2C A 



6P 



6/3 - 13C A 



1 + 



20 
13 



32 Cf_ 
13 

- 1 



Pa - 2C A 



ira s {m) < z 



,/3o-2C A 



0- 



(14) 



The {c} matching coefficients can be found or deduced from the results in Ref. ||. The 
{c?} matching coefficients are new. For the spin-dependent {d} matching coefficients there 
already exists an evaluation || but we differ with their result. This finishes the RG evaluation 
of the NRQCD Lagrangian at one-loop at 0(l/m 2 ). 

With the above results one can resum the large logs associated to the hard scale by 
running down the factorization scale u s up to the next relevant scale. 

Finally note that it is very important to know the basis of operators one has been working 
in NRQCD as well as in which gauge the calculation has been performed. In practice this 
means that one should make sure that cp and d vs have been computed in the same way in 
order to obtain the correct result. 



3 pNRQCD 

The above results are also a necessary step towards the RG improvement of pNRQCD when 
ma s S> Aqcd, which we consider in the following. 



6 



By integrating out some soft degrees in NRQCD one ends up in pNRQCD. This latter 
theory is defined by the cut-off v p nr = {v p , v us }, where v p is the cut-off of the relative three- 
momentum of the heavy quarks and is such that mv <C u p <C m and v us is the cut-off of 
the three-momentum of the gluons and light quarks with mv 2 <C v us <C mv. In principle, 
we do not rule out the option of correlating v v with v us in order to efficiently perform the 



renormalization group improvement at higher orders |TTJ. Nevertheless, at the order we are 
working with, we not need to specify any relation between u p and v us since the dependence on 
u p would be a sub leading effect. Therefore, in this paper, we will treat them as independent. 
The pNRQCD Lagrangian reads as follows: 

f / n 2 r> 4 \ 

£ p nrqcd = Tr St [ido - c k ^- + 04-^ - V}°\r) - '— + . . . S 

\ m Am 6 m m z I 



+0+ ^A,-c fc ^-yW(r) + ...^oJ 



+gV A (r)Tr {O f r • E S + S f r • E O} + #-y^Tr {O f r • E O + O f Or ■ e} 

1 

4 



— A Gl u G^\ (15) 



where we have explicitly written only the terms relevant to the analysis at the NNLL of the 
singlet sector; S and O are the singlet and octet field respectively. All the gauge fields in Eq. 



([Lq) are functions of the center-of-mass coordinate and the time t only. For a more extensive 



discussion we refer the reader to Refs. P, ITS] 



3.1 Potentials 

We now display the structure of the matching potentials VP , VP , VP and , which are 
the relevant ones to our analysis. 

1 ) Order l/m°. From dimensional analysis, VP(r) can only have the following structure: 

K (0) - -C f ^, (16) 
and similarly for the static octet potential: 

V* - (f - O,) (17) 

2) Order 1/m. From dimensional analysis and time reversal, VP can only have the 
following structure: 

VP C f C A D^ 



m 2mr 2 

3) Order 1/m 2 . At the accuracy we aim, VP has the structure 
VP C f D?l rl a1 , C,D®1 T2 , *C f D® 



18) 



m 2 2m 2 I r J 2m 2 r 3 m 2 

4ttC , / L'§ ) , iCfD i i ) s 1 C f D { sl 1 

+ 3m 2 S 6 W + -^m^^ • S + ~^7* Sl ^ (19) 
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where Si 2 (r) = 3r • <J\ r • <x 2 — <T\ ■ cr 2 and S = cr 1 /2 + cr 2 /2. 

The coefficients, ay s , D s , ... contain some lnr dependence once higher order corrections 
to their leading (non- vanishing) values are taken into account. In particular, we will have 
expressions like5 (3) (r)ln n r. This is not a well-defined distribution and should be understood 
as the Fourier transform of ln n \jk. Nevertheless, in order to use the same notation for all the 
matching coefficients, and since it will be sufficient for the purposes of this paper, resum the 
leading logs, we will use the expression 5^(r) ln n r, although it should always be understood 
in the sense given above. 



3.2 RG equations 

Formally, we can write the pNRQCD Lagrangian as an expansion in l/r(— 1/r, p) and 1/m 
in the following way: 

OO 1 OO 1 OO / 1 \ 

Vrqcd = E W + - E r»V™0™ + ± E r-V^O^ + O (^) , 

n= -1 m n=-2 m n=-3 Vm 7 

(20) 

where the above fields and parameters should be understood as bare. As for the renormalized 
quantities, we define V as the potentials and V as the (almost) dimensionless constants in it. 
The latter are in charge of absorbing the divergences of the effective field theory. Therefore, 
they will depend on v p and v us . Note that the dependence on v s of the NRQCD matching 
coefficients has to cancel in V since the new effective theory does not have any ultraviolet 
cutoff dependent on v s . This discussion completely fixes the procedure to obtain the RG 
equations of the potentials: by studying the UV behavior of pNRQCD it is possible to obtain 
the scale dependence of the potentials on v p and u us and the independence on v s trivially 
sets the v a scale (in-)dependence of the potentials. Being more specific, the potentials have 
the following structure: 

V(d(v p , v s , m), c(v a , m), u s , v us , r) = V(u p , m, v us , r) = V(u p , v us ) . 

This produces the following RG equations: 

v s -^-V(d(v p , is s , m), c(u s , m),v a , v us , r) (21) 



v — + v {—d\ — + v (—c\ — 

s dv s s \dv s I dd s \dv s J dc 



V(d,c,u s ,u us ,r) = 0. 



d ~ ( d \ d ~ 

v p —V{d{v p , v a , m), u s , v us , r) = u p I —d I g^ v (d, v a , v us , r) . (22) 

The first equation just reflects the independence of the potential on v s . At the practical 
level, with the accuracy we are working, it is equivalent to set v s = 1/r up to factors of 
order one. The second equation tells us that the dependence on v p is inherited from the 
(four-heavy fermion) NRQCD matching coefficients. 

One of our aims will be to obtain the heavy quarkonium spectrum with NNLL accuracy 
when Aq CD /(m« s ) 2 <^ ma 2 s . In that situation the leading order solution corresponds to a 
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Coulomb-type bound state and, leaving aside non-perturbative corrections, a perturbative 
expansion is licit. In order to achieve this goal we will need the RG improvement of the 
pNRQCD Lagrangian for the singlet sector with the same accuracy. Being more precise, we 
will need 

di/ dv du 

as well as (due to mixing) 

zyAy(o) _ a 2 z/-^-Va ~ a s , f t"Vb ~ a s , (24) 
dz/ dz/ df 

within an strict expansion in a s . 



Let us consider first the dependence on v p . Eq. ( P2"[ ) tells us that the dependence 
on z/ p appears due to the four-heavy fermion matching coefficients. These first appear at 
0(l/m 2 ). Therefore, we only need to obtain the equation v p -4j-V^ 2 ' ~ a 2 . In fact, at least 
at lowest non- vanishing order, only the delta potentials are dependent on u p . Of this type 
is precisely the leading dependence on v p of pNRQCD. It appears through (the iteration 
of) enough singular potentials when performing standard quantum mechanics perturbation 
theory. Explicit inspection shows that these kind of effects first appear at 0(ma®) in a 
perturbative computation of the mass (this argument is based on the singularity of the 
(iteration of the) potentials plus taking into account their leading non-vanishing power in 
the a s expansion). Therefore, we obtain 

^ (2) = ° + 0(4 (25) 

By using Eq. (^2|), this is equivalent to v v -^-d = + O(a^), which was needed previously in 
NRQCD. In principle, this result could also be proved by explicit inspection on the possible 
diagrams at the quark-gluon level that could give divergences proportional to v p . The final 
conclusion is that we can neglect any dependence on v p at the order we are working at in 
the potentials, i.e. V(is p , v ua ) ~ V{v us ). Therefore, we only have to compute the u us scale 
dependence. 



The z/^-scale dependence could be taken from the computation in [19], [18], [Hj, |3| (see also 
f2l"ll ) by keeping track of the dependence of the result on the different V. Let us note that we 
only need to do a one-loop computation in order to achieve the necessary accuracy (plus the 
already known two- loop singlet static potential). This should be compared with the one-, 
two- and three-loop calculations that seem to be necessary if the calculation is performed at 



the quark-gluon level as in Ref. |L2], [13], p!0| . 

Formally, the renormalization group equations of the renormalized matching coefficients 
due to the ^-dependence read 

Vus-^V = By{V). (26) 

From a practical point of view one can organize the RG equations within an expansion in 
1/m. 
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At 0(l/m°), the analysis corresponds to the study of the static limit of pNRQCD, which 
has already been carried out in Ref. [p0[| . We repeat the basic points here for ease of 
reference. Since V-\ 7^ 0, there are relevant operators (super-renormalizable terms) in the 
Lagrangian and the US RG equations lose the triangular structure that we enjoyed for the 
RG equations of v s . Still, if V-% <C 1, a perturbative calculation of the renormalization group 
equations can be achieved as a double expansion in V-i and Vq (for a similar discussion in 
the context of scalar A0 n -like theories see |22[|), where the latter corresponds to the marginal 
operators (renormalizable interactions)). At short distances (l/r> Aqcd); the static limit 
of pNRQCD lives in this situation. Specifically, we have V_\ = {ay s ,ay D }, that fulfills 
V-i ~ oc s (r) <C 1; Vq = a s {y us ) and V\ = {Va,Vb} ~ 1. Therefore, we can calculate the 
anomalous dimensions order by order in a s (v us ). In addition, we also have an expansion in 
V-\. Moreover, the specific form of the pNRQCD Lagrangian severely constrains the RG 
equations general structure. The result obtained in Ref. 



20 reads 



d 

V US -, OC Vs 

dv„ e 



dv us 



dv UR 



2a 



3 7T 

2a 



s v 2 



s \r2 



3 7T 



C f 



V 2 



Cf) a Vo + C f a v )j 



(Ca^ 



Cf a Vo + C f a Vs 



"us-. — V A 

du us 

d V 
Vus--, — Kb 

dv qis 



0. 



(27) 



At higher orders in 1/m, we only need to consider the singlet potentials. The same 
considerations than for the static limit apply here as far as the non-triangularity of the 
RG equations is concerned. At 0(l/m, 1/m 2 ), we have the following matching coefficients: 



V$ = pW, c k } and V% = {D\% D%, D%, D$., D^ s , and we obtain 



'(2) 



(2) n (2) n (2) n (2) 



>(2) 



,(2) 



dv,, 



d -Df] 



'dv u _ 

—D {2) 
dv,,* 



A 71 



16 a s 
~3~7 



a. 



V 2 r 2 
v A c k 



2„2 



3 7T 



'-Vic k 



% - C f ) a Vo + C f a Vs 



4 

c A 



2C f a Vs +(^-C f )a Vo 



Cf) av a , 



Cf ol Vo + Cfa Vs 



{21 



and zero for the other matching coefficients (in particular for the spin-dependent potentials). 
In a more formal way, Eq. (|28|) has the following structure: 



' dv. 



-V$ ~ VoV^V 2 , 



' dv. 



us 

d 



VoV^Tv^iV 2 . 



(29) 
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In fact, in general, we have the structure (Vffl = V m ) 
d 



v, 



'dv„ 



.y( n ) 



us 



Y v {ni) V {n2) ■■■V {n i\ with y 

/ j mi m.2 m j ' / j 

{rii}{mi} 



n, 



i=i 



3 

i=l 



(30) 



and one has to pick up the leading contributions from all the possible terms. 



Eqs. (pTD, flZoQ, ( Wn) an d (Pq) provide the complete set of RG equations at the desired 
order. By using Eqs. (j2lD and (|25|) , we obtain 



V = V{d{l/r), c(l/r), z/ s = 1/r, r) . 



(31) 



We now need the initial condition in order to solve the US RG equations, i.e. the matching 
conditions. We fix the initial point at v us = 1/r. In summary, we need to know the singlet 
static potential with O(a^) accuracy, the 1/m potential with O^af) accuracy, the 1/m 2 
potentials and the singlet octet potential with 0(a s ) accuracy and Va with 0(1) accuracy 
at v us = 1/r. They read 



av^r- 1 ) = a s {r- 1 )\l + (a 1 + 2 lE (3 ) 



47T 



+ 



7T 



lE (4aiA) + 2/3 x ) + — + He] Po + «2 



16 7T 2 















r 




= a s 






y 




= a s 










= a s 


r- 1 )(2 + c D (r- l )-2c F (r- 1 )) 



+- 



7T 





= a s 


[r 


l ) 






[r~ 


- l ) 




3 






= a s 


[r~ 


l ) 


ttVo(^ _1 ) 


= a s 


J~ 






= 1, 







d vs {r- 1 ) + 3d vv (r~ 1 ) + —{dUr^ 1 ) + 3d sv (r~ 1 )) 
(d«,(r _1 ) + C / d BW (r~ 1 )), 



* v ' 2nC f 
(c s (r- l ) + 2c F (r- 1 )) 

,2 / -1\ 



(32) 

where (3\ = 34/3C A — 4CfT F nf— 20/3C^T F n/ and the values of a x and a 2 have been computed 
in Ref. [f23fl . We now have all the necessary ingredients to solve the RG equations. 

Eqs. fl2"?D and ( p5|) give rise to sub leading effects within an strict expansion in a s . 
Therefore, we can approximate them to (if not displayed the RG equation remains equal) 



d 



' dv q . 



-a Va 



2 a s (u us ) 
3 



7T 
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US 7 ^Vo 

du us 



'us 

ly us 



J_ D &) - 8 a s (u us ) C A ! 
us dv us Dl > s - 3 vr 2 °' (r j ' 

d ( 2 ) 16a s (z/ us ) /C A \ u 



D% = ^-^l l^-C f WO- (33) 



'di/ us ,s 3 vr V 2 
We can finally obtain the RG improved potentials for the singlet: 

<*.(*.) = «K.(r-) + f^(r- )bg (|g), 

= D ( 1)(r - 1) + i|(^ + C/ ) a , V1)log (^gl- 

Cg(^) = <2r 1 ) + ^ia,(r-')log(^), 

3po V^sl^s)/ 



6p \ A J \ OL s \V us ) ) 

D%(u ua ) = flg^r- 1 ), 

= flSUr -1 )- (34) 

This completes the RG evaluation of the pNRQCD Lagrangian at NNLL (as far as the singlet 
is concerned). 



4 Heavy quarkonium spectrum 

In the situation Aq CD /(ma s ) 2 <C ma^, the heavy quarkonium behaves as a Coulomb-type 
bound state and the perturbative corrections can be computed in a systematic way. From 
the potential-like terms, we obtain the following correction to the NNLO energy expression 
(the derivation of this result would go along similar lines to those in Ref . |TS[ ) : 



, C f 5l ° 



E n a s 



3n 



16 

To 



4Cpi 8 s i 
3 n 



2C A 

_Ca 
2 



24) 



Q+4C A C f , } x 
2 / n(2i + l) 



+ (^TTT) 



1 



n- 



log 

3 

27ra„ 



ct s {y us ) 



0:.s 



-2C A 



C A 



z -/3o _ Z -2C A 



log 



(35) 
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where E n = —mCja 2 /(4n 2 ), the scale v s in z and in the NRQCD matching coefficients has 



J f 

been fixed to the soft scale v s = 2a" 1 where 2a" 1 
soft scale v s = 2a" 1 unless the scale is specified, and 



mCf a s (2o n 1 ) 



. a, is also understood at the 



(1 + 1) 



2/ - 
-4 (z 

27T3{ 4 ' 2, + 3 H 



!) + (, 



{ {4(2/ - 1) (z 

l) + (V 2Ca - 1 



-2C A 



0} 



-2C A 



0} 



(36) 



Z + l. 



Eq. fl3"5| ) gives all the 0(maj (a s ln) n ) terms for n > 1 of the heavy quarkonium mass, where 
In stands either for ln(a s ), arising from the hard scale, or for Yti{ma s / v us ), arising from the 
ultrasoft scale. After adding to Eq. fl35|) the NNLO result with the normalization point 
at the same soft scale, v s = 2a" 1 , that we have used here, the complete (perturbative) NNLL 
heavy quarkonium mass is obtained (note that for the LO result the three-loop running a s 
has to be used). 

The v us dependence of Eq. fl3"5p cancels against contributions from US energies. At the 
next-to-leading order in the multipole expansion, the contribution from these scales reads 



5 VS E 



3N, 



■Tp / ^(n,/|re^"- Ho) r|n,/)(E a (t)0(t,O)^ J E b (O))(^ 
Jo 



(37) 



P 

where H Q = Ck h and v us is the UV cut-off of pNRQCD. Then, the total correction 

m 

to the energy reads 

8E nM = 5 pot E n ^(v us ) + 5 vs E n ^ us ) . (38) 

Different possibilities appear depending on the relative size of Aqcd with respect to the 
US scale ma 2 s . If we consider that Aqcd ~ ma 2 s) the gluonic correlator in Eq. (|37|) cannot be 
computed using perturbation theory. Therefore, in a model independent approach, one can 
leave it as a free parameter and fix it with experiment at some scale v us (since the running 
of Eq. ( |3"7| ) with v us is known, one can then obtain its value at another scale). 

If we consider that ma 2 s 3> Aq CD , Eq. (|37| ) can be computed perturbatively. Since ma 2 is 
the next relevant scale, the effective role of Eq. ( ]3"7| ) will be to replace u us by met 2 (up to finite 
pieces that we are systematically neglecting) in Eq. (p5]). Then Eq. (|38| ) reduces to Eq. ( p5[ ) 
with v us ~ ma 2 . In particular, we take v us = —E n . As expected, Eq. ( |35|) with v us = —E n 
reproduces the already known 0(ma 5 s In a s ) correction ||18| (see also 



situation one is assuming that 



A 



3 0) 



QCD 



mat 



< 1, 



Since in this 
one can expand on this parameter. Therefore, 



non-perturbative corrections can be parameterized by local condensates. The leading and 
next-to-leading non-perturbative corrections have been computed in the literature |26|, |27l . 



There already exists an evaluation |10[ within the vNRQCD framework |TT| , 12 , |l"3f of 
the RG improved Heavy Quarkonium mass when Aqcd "C ma 2 . We agree for the spin- 
dependent terms (since we agree with the spin-dependent potentials computed in Ref. ||12|| , 
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see Eq. (|34"D) but differ for the spin- independent ones. We note that the disagreement still 
holds if we consider QED with light fermions (Cf — > 1, Ca — > 0, Tp — » 1). Agreement is 
found for QED without light fermions (Cy — > 1, — ► 0, n/ — > 0, Tp — > 1). 
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